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I. INTRODUCTION 



What is now referred to as the P-L-K (Poincare-Lighthill-Kuo) method 
or the Lighthill method of strained coordinates was first proposed by 
M. J. Lighthill in 1949 [1] as a perturbation method for obtaining uni- 
formly valid approximate solutions for certain classes of ordinary and 
partial differential equations. It was then used successfully by many 
authors (see Lighthill's 1961 lecture [2] fora sample list) although there 
was a growing feeling of uneasiness on the part of some workers because 
of cases of partial differential equation problems where the method gave 
some wrong answers, even though the method looked right (see Van Dyke 
[ 3 ] , Levy [4 ] and Lin [5 ] ) . 

Because partial differential equations are too hard, mathematical in- 
vestigations have centered on the ordinary differential equauon problems, 
and in particular, the model equation originally used by Lighthill. 

(x + ey) q (x)y = r(x) 0 <; x £ 1 (1) 

y (l ) = b 

Following a conversation at a cocktail party Wasow [6] proved that 
Lighthill's method worked for (1) with an added hypothesis. With Wasow' s 
paper the purely mathematical interest died, except for isolated papers 
which went basically unnoticed. The essential ones are Temple's lec- 
ture at the 1958 International Congress [7] where he proposed a modi- 
fication which amounts to a different motivation of a similar method, 
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Pritulo's paper [8] which reduces a portion of the problem to an algebraic 
equation as opposed to a differential equation, and Takahasi's paper 
[9] which generalized (1) to include higher order non-linearities and 
proved the convergence of Lighthill's procedure in the manner of Wasow. 

Since 1968 however, there has been considerable activity in study- 
ing the validity of applying Lighthill's method to the ordinary differential 
equation case. The activity started with Comstock's paper [10], based 
on a conversation (at a cocktail party) with Lin, which showed by a 
series of examples that Lighthill's method was quantatively inaccurate 
in the more general case, 

(x n + ey) ^ + q (x)y = r(x), n > 1 (2) 

and, in addition, the method was also qualitatively inaccurate if Wasow' s 
extra criterion did not hold, although the method appeared to work in all 
cases. Since that paper, several papers have appeared trying to patch 
up these difficulties. Usher [11] rediscovered a portion of Pritulo's 
result [8] and showed where Wasow' s criterion seems quite reasonable. 
Burnside [12], starting from Temple's approach, showed that an initial 
change of independent variables to reduce the exponent n in (2) to one 
would solve some of the inaccuracies pointed out by Comstock [ 10]. 
Independently Melka [13] arrived at a similar, and slightly more general, 
result. Also, in an as yet unpublished work [14], Melka suggests a 
different and potentially interesting approach. And most recently Su and 
Liu [15] finally attempt to analyze the problem in the framework of a 
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general asymptotic expansion. 

The purpose of this paper is to combine these results to illustrate 
the different approaches and their inherent problems, correct some errors 
and to suggest some new results. 

II. THE P-L-K METHOD AND ITS RESTRICTIONS 

The original model of Lighthill [1] was the problem 

( x + ey) ^ + q(x)y = r(x) 0 <; x ^ 1 (1) 

y(l) = b > 1 (2) 

q(o) t 0 (3) 

This last condition is one to which we will return later. If one tries to 
find y as a series in e 

2 

y ~ y (x) + cy (x) + e • y (x) + ... 
o 1 2 

then the equation for y Q is 

dy o 

X + q(x)y Q = r(x) (4) 

whose homogeneous part has a regular singular point at the origin. The 
original equation, as a simple phase plane analysis shows, does not 
have a singularity at the origin, and thus the perturbation series has a 
false singularity. (It is the condition (3) which makes this singularity 
a regular singular point.) The false singularity is disturbing enough, but 
in addition it is easy to see that the solutions for y^x), y 2 (x), etc. are 
even more singular, due to the term of y dy/dx. 
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Lighthill's idea was to move the singularity back out of the domain 
of interest by introducing a new, slightly stretched, independent vari- 
able z by the implicit equation. 



x = z + 



ex^(z) 



2 

€ X r 



(z) + 



(5) 



and then look for a solution y in the form 

2 

y = y (z) + e y, (z) + e y 0 (z) + ... (6) 

o i l 

Utilizing (5) and (6) in (4) and expanding the functions r(x) and q(x) 

about the point z, one generates a single sequence of differential equa- 
th. 

tions. The i equation is a single linear first order equation in both 

dy./dz and dx./dz, in terms of the lower order y. and x. . Clearly one 
ii ii 

must generate a second sequence of equations to determine both the y^ 

and the x^. Lighthill's choice is best phrased by Van Dyke [3] vvho said 

"higher approximations shall be no more singular than the first. " In this 

th 

instance one groups all the non-homogeneous terms in the i equation 
which could contribute to making the solution for y^(z) more singular than 
y i _^(z) into a group including all the x^z) terms, and set this equal to 
zero, thus creating a differential equation forx^(z). This necessitates 
chosing a boundary condition for x^(z) as well. One such choice is to 
make x = 1 correspond to z = 1 , i . e . , 

x.(l) = 0, (7) 

although this may not be the best choice. 

We reemphasize the choice of grouping of terms to create a 
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differential equation for x., . One only need put the singular terms into 

this equation. One can go so far as to put all the non-homogeneous 

terms into the equation for x (z) . Then, with the boundary conditions 

y.(l) = 0, i the linear homogeneous equations for y will make all 

the y . s 0 and equation (6) converges very nicely, consisting only of 

y (z). 
o 

Wasow in a paper [6], which was extended and corrected 11 years 
later by Sibuya and Takahasi [16], showed that the above procedure, 
choosing y . - 0 i ^ 1 , led to a convergent series expansion for the 
x,,(z), provided q(0) ^0 and also 

y Q (z) q(z) - r(z) 4 0 (8) 

0 < z £ 1 . 

This inequality has since been referred to as the Wasow criterion, and 
did not appear in the original Lighthill work. The paper of Sibuya and 
Takahasi went a bit further and asked whether this convergent series 
for x-z(z, e) could be inverted for z. They prove that in addition if 
there exists any analytic particular solution cd(z) to the zero order equa- 
tion with the property od(l) ^ b, then the series is invertable. If 
W(l) <b, the point x = 0 corresponds to a positive value of z and there 
is no difficulty. If to (1) > b then the solution to (1) has a genuine singu- 
lar point in the domain of interest and the point x = 0 corresponds to a 
negative value of z . In a subsequent paper Takahasi [9] extended these 
results to the equation 
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( 9 ) 



a , c 

{x + L t™ P m (x,u)} ~^ + q(x)u = r(x) + L € R m (x,u) 

where the P m and the R m are polynomials. The theorems are essentially 
the same as the previous paper by this author [16]. 

Thus the work of Wasow and his successors have led to the following 
result: 



For the equation 



du 



(x+ e P(u)) — + q(x)u = r(x) 
u(l) = b 



( 10 ) 



Lighthill's procedure of expanding both u and x in terms of a slightly 
strained variable z, given by 



n , . 

x = z + S e x (z) , 
1 



(ID 



subjecting the x^(z) to the differential equation naturally generated by 
using series expansions in (10) and grouping the singular terms , will 
work provided certain conditions hold. They are 



i) q (0) 4 0 



ii) u q (z) q(z) “ r ( z ) ^ 0 in 0 < z ^ 1 



( 12 ) 



iii) b is large enough 



Three questions immediately arise. Can the procedure be extended 
to more general equations than (10)? Secondly, must the straining be of the 
form (11) or would a more general straining work? And lastly, what is the 
significance of the restrictions (12) listed just above? 
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Temple's approach to Lighthill’s problem [7] certainly suggests 
that a much more general problem could be attacked. Temple framed the 



problem in the following way. Rewriting (10) as 



du _ q(x)u + r(x) 
dx x + eu 



(10a) 



one can view this as the result of combining two ordinary differential 



equations 



du 



d4 

dx 



= q(x)u + r(x) 



= x + eu 



(13) 



Now the differential equation for x is chosen ab initio. Because of the 
analyticity of the right hand sides of (13) in x, u and e, a series ex- 
pansion for x and u in powers of e is assured to converge The variable 
£ in (12) is related to z in (11) by £ = log z. 

There is a slightly more natural way to get Lighthill's stretched 
variable z by writing (10a) as 

z^T = q(x) u + r(x) 

dz (13a) 

dx 

z— = x + e u 
dz 

However these equations have an obvious singularity, where as (13) 
does not. There is now an elementary jump to the equation 

[x k + eg(x,u,e)] h(x,u,c) = 0 (14) 

which can be written 
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-htx.u.c) 

dx k , 

x + e g(x,u,e) 



( 15 ) 



and Bellman [17] makes this jump, saying that Temple's method would 
work just as well here. As we shall see, this conjecture is false. 

Temple's approach actually presents several difficulties, but it also 
suggests a couple of positive conjectures. The first difficulty is that the 
differential equation for x is selected with no reference to deleting any 
singularities for the expansion of u. This removes the original motivation 
for the problem. In contrast to the series used by Wasow, both the 
u series and the x series generated by this split will, in general, be 
infinite series. The flexibility to make the u expansion a finite number 
of terms is gone. A second difficulty is that the Wasow restricts s 
(12) are all on the equation for u, which now does not appear to have any 
problem. The biggest difficulty, however, is that the singularity involved 
in the equation is hidden. The equations (13) are actually analytic in 



u, x and e , provided 2 is finite . But the zero order solution for x is 

z 



x = e 
o 



(16) 



and the compact domain 0 ^ x ^ 1 is mapped into the semi-infinite domain 
-® ^ z <0, and unfortunately it is the point at infinity which is the point 
of primary interest for most of the problems. We will return to this point 
later. 

Another splitting of the equations, starting from the Lighthill approach. 
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was devised by Pritulo [8] and later again by Usher [11]. A simple 
observation reduces the problem of finding x. (z) to an algebraic problem 
rather than a differential equation problem. This observation is the 
following theorem [8]. 

Theorem: If an ordinary perturbation expansion of u alone has a (non- 
uniformly convergent) expansion 



U ~n=0 f% n ,xK 

then the Lighthill expansion for u 



is given by the formula 



(17) 



where 



with 



where 



u ~ L e 11 U (z) 
n=0 n 



n-1 



U n (z) - u„(z) + E C n (z U u (z) 
n=0 



c (z) = x (z) 
o 1 



c p ( z) = k = (2k - p) x k + i (2) c P -k (z) 



n 



x = z + £ e x. (z) 
1 i 



(18) 



(19) 



( 20 ) 

( 21 ) 



The power of this theorem is that these equations (19-21) hold for 
any choices of the x^(z). Thus we can make certain choices. By hypo- 
thesis the u (x) form an increasingly singular set, so the U (z) do also, 
n n 

However, in (19), the x^(z) can be chosen successively so that they 
subtract off the increasingly singular part, so that the U n (z) can be made 
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uniformly convergent! 

Prit ulo extended his formula to higher order differential equations, 
but he made no applications of his results. He was also apparently not 
aware of Wa sow's paper. 

Several years later Usher made a similar observation [11]. Usher 
started from the equation 

“j^= f(x,y,e) (22) 

where y may be a vector and assumed that (22) was analytic in e so that 
it may be expanded in the form 

^•=f (x,y) + f f (x y)+ ... (22a) 

dx o 11 



Looking at the model equation (1) in this way 



dy _ q(x)u + r(x) 
dx x + e u 



0 sx £ 1 



1 



= “ [q(x)u + r(x) ] + 77 [q(x)u + r(x)] + . . 

X X 



we see that the expansion is invalid in the given domain. However, 
Usher only wanted the first order terms in u and x, and the lack of 
analyticity in e does not affect these calculations. Usher then looked at 
the zero and first order terms, in light of Wasow’s results. We write 
them out, in Pritulo's notation 



U (z) = u (z) 
o o 

du 

Uj (z) = Uj (z) + Xj (z) 



(23) 
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But since, for the model equation, 



du 

dz 



= - [q(^)u o 



+ 



r(z) } 



we have 

X 1 

U (z) = u (z) + — [q(z)u + r(z)]. 
1 1 Z o 



( 24 ) 



Separating u^(z) into its very singular and mildly singular parts 
Uj(z) = Uj m (z) + u j s ( z ) we see that there is a solution of (24) for x^z) 
which leaves U^(z) only mildly singular provided that Wasow's criterion 
holds. Thus it was Usher who realized that the restriction ii) in (12) 
(Wasow's criterion) is necessary criterion for a finite x^(z) throughout 
the range 0 <z s 1. And equation (19) shows that z ^(zju^ + r(z)] 
is the coefficient of each succeeding xdz), so this restriction applied to 
each succeeding term. 

Does the singularity in the expression forx^ invalidate the use of the 
method? Apparently not , to first order, as the following example shows. 
We consider a problem which is explicitly solvable, namely 

(x + eu)^- =-u+2bx Osx^l 

dx 

(25) 



u(l) = B = b(l f 77 ) 



By simple integration, since the equation is exact, we have 



with 



1 2 2 

u = e { -x + Jx + 2 e (c + bx )} 

- k b 2 , 

c = hr} + e 2 d + V> 



(26) 
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Solving this problem using Usher's approach, we write 

u(x) = u (x) + eu. (x) + ... 
o 1 

u (1) = B (27) 

o 

u^l) = 0 i ^ 1 

this gives 

u (x) = bx + b«x (28) 

o 

2 

u^x) = j {(1 + rf) x 1 -rj 2 x 3 -x} (29) 



and the Wasow criterion term is 

-u + 2bx = b(x - ^) (30) 

o x 

We recognize that u n (x) is more singular than u (x) , and also that (30) 

1 o 

vanishes in the domain 0<x<lif0<?7<l. Using Usher's approach we 
rewrite 



U (z) = b(z + n ) 
o z 

u l(z) = j {(1 + rj 2 )z 1 -z -r/ 2 z 3 ) 

= U x (z) - Xj (z)b { 1 - n 2 } 

z 

There is some flexibility in the grouping of (29a) to create U^(z). 
choose to make U^(l) = 0 and make the choice 

U x (z) = | (1 + r7 2 )(z" 1 -1) 



(28a) 

(29a) 



We 



(31) 



so that 

x l(z) = j f r? 2 ( z 3 -1) + (z-l)){l - } 1 (32) 

z 

and x^ (z) is indeed singular at z = Jrf. 
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However, what information do we wish to know? We claim that at 



this stage the only important questions are a) does the qraph of U 

o 

approximate that of the true u, and b) what is the value of u(0) predicted 

by our approximation. Ignoring the singularity in x^(z) we can graph the 

sum U (z) + c U n (z) , and since u is no more singular than U (z) it will 
o 1 o o 

not alter the graph. Now U q (z) is strictly positive and has a minimum at 
z = Jr) where its value is U q C Jr\) = 2b Jr\. But an investigation of the 
exact solution (26) shows that u(x) is strictly positive, and has a zero 
slope at x q = Jr\{\ + eb) 1 , where its value is u(x q ) = VoJr\ (1 + eb) *. 
And if we assume that x = 0 still corresponds to z very small then we 
can solve the equation , 

0 = x = z + ex^z) = z + ^ {rj 2 (z 3 -l) + (z-1)} {l - n 2 } 1 (33) 

z 

for z, at least approximately. We obtain 
z(0) b|- 

so that 

u (x = 0) „ (34) 

o z e 

as compared with the exact value from (26) 

u(0) = J2f , c = b„ + £§- (1 + 17 2 ) (35) 

Thus these features of the approximation fit quite well. 

Is our assumption that z is small corresponds to x = 0 valid, in 
view of the singularity in x^(z)? A study of x^(z), equation (32), shows 
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that this is valid. An elementary analysis of the numerator of (32) shows 

that x^(z) has two positive zeros, one at z = 1 by construction, and 

one at some smaller value of z, depending upon r). The details of the 

behavior depend upon whether this zero is greater than or less than the 

singular point z = >Jr\. We sketch the graph for the two cases. The 

critical value where the zero and the singular point coincide is given 

dx 

by looking for the switch or sign of — near Jr\ which occurs at 

1 + v 2 = 2fn (36) 

or 

rj - .296 (36a) 
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The fact that this example gives a number of right answers despite 
the extreme distortion of the "stretched" coordinate encourages us to 
ask whether the distortion can be corrected. To see that it can, we 
return to a comment made earlier concerning the flexibility on choosing 
x^(z) and U^(z). In particular, there is no necessity for making x and 
z coincide at x = 1. (Doing so vastly simplifies the assigning of the 
condition u(x = 1) = B, of course.) In particular one can make x and z 
coincide at the singular point x = Jy\, as we show for this example. 

Going back to (29a) we write 
2 

Uj(z) = — p{(l + ? 7 2 )z 1 + 2z - 3z - rj 2 z 3 + Vt) - 4 Vrj) 

2 2 __ < 29b) 

= — {(L + 77 2 )z _ 1 + 2z - Vt?} - Y {t? 2 z 3 + 3z - Vt?} 

Now we have 

Xj(z) = j {t) 2 z 3 + 3z -Vrj} {1 - 7?z 2 3 (37) 

which has the limit 

x 1 ( v /^) = 0 (38) 

Obviously all higher order x^(z) can be treated in a similar fashion. And 
we observe that we have made the zero of x^ a simple zero at z = Jr\, so 
x^ is negative to the left of x = ^/r^and positive to the right, as it should 
be . 

To obtain the correct boundary conditions for u in the series form 
(27) we can not take the u^(z) = 0 at x = 1 , as we did. This means 
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that we must go back one step further and take the indefinite integral 

2 

for (29b), replacing (1 + t) ) by c^. Then we have 

h 2 -1 

U 1 (Z)= 2 {C 1 Z +2z -4^} (39) 

To determine c^ we must examine the equation 

1 = z + exj (z) (40) 

and solve for z. We get 

z(l) = 1 + edj + 0(e 2 ) 

Then we look at 

+ °(€ 2 ) = B 

= b[l + ed^ + t)( 1 + ed^) (41) 

b 2 -l 

+ €~2 [Cjd + edj) + 2(1 + edj) - 4 „frf\ 

Equating powers of e we get the condition for c^ to be 

Cj = Afrj- 2 -^d^l-rj) (42) 

Higher order terms may be computed i n an analogous fashion. 

We now , conjecture the theorem: 

Theorem: In the case that the zero order term in the Lighthill expansion 

for (1) has a single zero slope at some x in the interior 0 < x < 1 , then 

o o 

the series expansions may still be made to converge by making x = z at 
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X Q . The theorem would be much more difficult to prove than the Wasow 
result, because by necessity both the x series and the u series are 
infinite series. 

The fact that it might be advantageous not to make z = x at x = 1 
was observed by Melka [13] in his thesis, but the reasons for making 
any particular choice were not discussed in any detail. 

We now turn to the question of the restriction q(0) ^ 0 (Eq (3)) in 
Lighthill's original problem. Wasow also makes the restriction in his 
theorem (12i). Lighthill's choice is based upon the character of the 
singular point in the zeroth order equation near x = 0, i.e. , 



the case q (0) >0 is of different character than q(0) < 0. However, 
several recent papers have, without saying so explicitly suggested that 
the restriction is more fundamental than Lighthill suggested. 
Comstock's [10] examples were of the form 



For n > 1 the zeroth order equation has the same character as Lighthill's 
case for q > 0. However, even in the cases where Wasow’ s criterion 
[ 1 2 ii ] held and the P-L-K expansion of (43) had the correct qualitative 
behavior, the P-L-K solution to (43) give 



x — + q(0)u = r(0) 



has a solution which behaves like u « A q x m near the origin, and 



r n du , n-1 

[ x + eu] — + nx 



u = mx 



(43) 



u(l) = b > 1 
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as opposed to the exact expansion 



u(oj - ySS) 

€ 

The error is a multiplicative factor of Jn . 

Carrier has given a famous example [18] to show that the zero 
order equation need not be linear in x. He studied several equations of 
the form 

(x 2 + cu)^ + “U = r(x) (44) 

ax ol 

For the case r = 0 this equation is exactly integrable in terms of Bessel 
functions and one obtains the implicit solution 

i [J, (a/eu) + A y. (a/eu)} 

x = — (cu) 2 (45) 

[j (ooJTu) + Ay (a/Ta) 

o o 

with the asymptotic approximation 

u(0) « [c j] _1 (46) 

The P-L-K method works quite satisfactorily to give this result also. 

Burnside [12] in trying to correct the discrepancy of Comstock's 

examples for n 1, took the approach that the essential feature is to 

du 

linearize the x dependence in the coefficient of — . He always makes 

dx 

the change of variables 

z = x° (47) 



so that (43) becomes 



, > du m 

(z+ €U dz + u = n z 



-1 



(48) 
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and then he expands z in a Lighthill expansion. In this way he con- 
cludes that the Jn discrepancy comes from the neglect of n terms of the 



f Vi 

same order due to an expansion of a term raised to the n n power. We 
note that (48) has q(0) ^0 whereas (43) does not 

Melka [13], independently of Burnside, tried to eliminate the dis- 
crepancy in Comstock's examples by another change of independent 
variables. He made the choice for a Lighthill expansion of 
k k 00 



x = z + S € x (z) 
1 n 



(49) 



where k is to be chosen in some fashion which is never specified except 
by example. For the Comstock examples he noted that a boundary layer 
approximation near x = 0 suggested that k = n is the "logical" choice. 
This choice will also eliminate the discrepancy. 

Melka also considered a generalization of the Carrier problem 

(x 11 + eu)^- + ax P u = 0 (50) 

P_|~ ^ 

in the case n = 2 (p + 1). This equation is a Ricatti equation for x 
and then Melka claims that an obvious choice in (49) is k = p + 1. 

(Note for Carrier's example this gives just the ordinary Lighthill 

expansion.) His analysis again is based on knowing the exact solution. 

We suggest that there is another basis for a choice of change of 

coordinates , a choice that will make the Lighthill method valid. Our 

choice is the following: given the equation 

(x 11 + eu) “ + q(x)u = r(x) (51) 
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with q(x) ~ ax P as x -* 0, first make the change of variables 

P+1 

y = x 

so as to make q(y)-» a, a constant ^0 as y -» 0. In this new coordinate 
system Lighthill's restriction (3) is met. In all the known exactly in- 
tegrate cases the author can find, the subsequent Lighthill approxi- 
mation is valid. In other cases the Lighthill procedure appears to work, 
but we are unable to check our answers against a known solution. It 
is the property q(0) ^ 0 that we feel is the important criterion in the 
choice of variables, not the resultant linearity or non-linearity of the 
coefficients . 



III. GENERALIZATIONS AND CONCLUSIONS 

The original interest in the PLK method was for applications to 
partial differential equations. There are essentially no theoretical re- 
sults concerning the application to partial differential equations. The 
primary results are examples and counterexamples. Lin [5] and Van 
Dyke [3] have concluded that this method works quite well in hyperbolic 
equations, where the characteristics become the "strained" coordinates. 
For many years it was recognized that the application of this method to 
elliptic equations usually led to erroneous results [2], [3]. However, 
in recent years two papers have appeared deriving successful approxi- 
mations for elliptic problems where previous workers failed [19], [20]. 
We will sketch the results obtained. The problem is that of flow around 
a thin airfoil, which can be described by finding a harmonic function 
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<p(x,y) satisfying the following conditions: 



< p -* Ux + 0(1) 



2 2 

as x + y -♦ co 




on the surface y = f(x) 



x 



2 — 

where f(x) = ch(x). In particular one choice for f (x) is f(x) = + c (1 - x ) 2 , 



a very narrow ellipse. The structure of the mathematical problem is 
significantly different from the ordinary differential equation problems, 
but the same heuristic is still present. For e ^ 0 the problem has no 
singularity in the domain of definition (the singularity is the "source" 
and "sink" inside the thin wing which "create" the wing). In the limit 
e = 0 the singularity is at the edge of the domain. An attempted per- 
turbation expansion for the solution in powers of e times a function of 
the complex variable z exhibits the familiar problem of increasingly 
singular terms. 

What is the straining of coordinates which is appropriate to this 
problem? After several erroneous attempts by several people, Hoogstraten 
[19] recognized that a key was the preservation of analyticity. He asks 
that 



map the body exactly onto the -1 £ 1. This forms his set of 

"boundary" conditions for the z , in a fashion very analogous to setting 

n 

z = 1 <=> x = 1 in the ordinary differential equation case. 

Martin [20] applies a somewhat different criterion, again based on 
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analyticity. He considers the same problem but looks at the expansion 

V 

CO f 

z = r? + L € Z (fi 
n=l 

as a case where one can use the Lagrange inversion formula. He then 

adopts Pritulo's formula for the expansion of the dependent variable. 

Thus he is not thinking in terms of a choice of boundary conditions for 

the unknown stretching, but of an algebraic criterion for this choice. 

He makes the following choice. Looking at the sequence of equations 

for the dependent variables choose ^(17) = where is the lowest 

order solution, and the a are (unknown) constants. It is easy to show, 

n 

for linear equations, that this can always be done. The result is that the 
solution for <p sums exactly to 

cp(r 7) = g(e) </? 1 (77) 

where g(c) = La c n . This gives a system of algebraic equations for 
n 

the z (77) , as in Usher's [ 11 ] paper. Martin gives no criterion for the 
n ' 

choice of the constants a^, and he gives examples where some choices 
give the right answers and other equally reasonable choices give non- 
uniformly valid answers. It is not at all apparent how to chose correctly 
since his choices are based on algebraic convenience and knowing the 
answer rather than some physically justifiable a priori criterion. The 
examples give no firm clue as to a choice. 

We may then conclude the following . For the ordinary differential 
equation problem, the hypotheses of the original Wasow proof are not 
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only sufficient but probably also necessary, with one exception. It 
appears that we may be able to correct for one zero in the function 
y^(z) q(z) -r(z) (see (8)) by a change in the stretching. Thus Usher's 
"necessary condition" are not as necessary as he thought [13] . For 
the generalized equation (2) it appears that the Lighthill procedure will 
work with a somewhat different change of variables, namely (49) where 
k is chosen to ensure that the criterion q(z) ^ 0 when z-* 0. For the 
partial differtial cases, in particular elliptic problems, it appears that 
the essential feature, for linear problems, is to make an analytic change 
of variables, although there is not yet any analogue of the Wasow paper. 
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